Effective Field Theory for Low-Energy Two-Nucleon Systems 
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We illustrate how effective field theories work in nuclear physics by using an effective Lagrangian in 
which all other degrees of freedom than the nucleonic one have been integrated out to calculate the 
low-energy properties of two-nucleon systems, viz, the deuteron properties, the np ^Sq scattering 
amplitude and the Ml transition amplitude entering into the radiative np capture process. Exploit- 
ing a finite cut-off regularization procedure, we find all the two-nucleon low-energy properties to 
be accurately described with little cut-off dependence, in consistency with the general philosophy of 
effective field theories. 
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Effective field theories (EFTs) have long proven to be 
a powerful tool in particle and condensed matter physics 
[Q,^, so it is quite natural that a considerable attention 
is nowadays paid to the role of EFTs in nuclear physics 
where phenomenological approaches have traditionally 
been tremendously successful. Some authors have fo- 
cused on nucleon-nucleon interactions and two-nucleon 
systems ^-|^ while some |^ on many-body systems in- 
cluding dense matter relevant to relativistic heavy-ion 
processes and compact stars. One of the most spectac- 
ular cases was the recent chiral perturbation calculation 
of the radiative np capture at thermal energy jij with an 
agreement with experiment within 1%. What was cal- 
culated in Q was however the meson-exchange current 
corrections relative to the single-particle Ml matrix ele- 
ment with the latter borrowed from the accurate Argonne 
Wis ]l0C phenomenological two-nucleon wave function. In 
this respect, one cannot say that it was a complete calcu- 
lation in the framework of the given EFT, namely, chiral 
perturbation theory (ChPT) although it was following 
the strategy of Q of using ChPT for computing "irre- 
ducible graphs" only. 

The purpose of this Letter is to supply the "missing 
link" that can render Ref. jij a "first-principle calcula- 
tion," that is to obtain the single-particle Ml matrix el- 
ement within the framework of EFTs In so doing 
we will compute the static properties of the bound np 
state (the deuteron) and the np scattering amplitude in 
the "^Sq channel. The results come out to be in a sur- 
prisingly good agreement with the data, offering a first 
glimpse of how EFTs work in nuclei. 

Since we shall be interested in very low-energy pro- 
cesses with the energy scale E <C ~ 140 MeV, we 
will integrate out all massive fields as well as the pion 
field j|] , leaving only the nucleon matter field which can 
be treated in heavy- fermion formalism (HFF). Since the 
anti- nucleon field also is integrated out in HFF, there 
are no "irreducible" loops (there will be, however, "re- 
ducible" loops to all orders in solving Lippman-Schwinger 
equation) and the EFT becomes non-relativistic quantum 



mechanics where all the interactions appear in the poten- 
tial. Now the np states we shall study are all very close 
to the threshold: they are either weakly bound i^Si) 
or almost bound (^^o). Bound states are not accessi- 
ble by perturbation expansion and the scattering state 
with a large scattering length a has a small scale a~^, 
making the convergence of EFTs highly non-trivial. We 
shall circumvent these difhculties by summing "reducible 
diagrams" - which amounts to solving Schrodinger (or 
Lippman-Schwinger) equation and using a cut-off regu- 
larization instead of the usual dimensional regularization. 

Due to the nonperturbative nature of the Schrodinger 
equation, unlike perturbative cases, it does matter which 
regularization scheme one uses in effective theories. Ka- 
plan, Savage and Wise and Luke and Manohar 
have found that with the dimensional regularization, the 

EFT breaks down at a very small scale, Pcrit — \J fo'^ 
large scattering length a, where r^ is the effective range 
and that this problem cannot be ameliorated by intro- 
ducing the pionic degree of freedom. As pointed out by 
Beane et al |^ and Lepage [|2j, the problem can how- 
ever be resolved if one uses a cut-off regularization. In 
effective theories, the cut-off has a physical meaning and 
hence it should not be taken to infinity as one does in 
renormalizable theories . In fact the strategy of effec- 
tive field theories is such that one should not pick either 
too low a cut-off or too high a cut-off: if one chooses 
too low a cut-off, one risks the danger of throwing away 
relevant degrees of freedom - and hence correct physics 
- while if one chooses too high a cut-off, one introduces 
irrelevant degrees of freedom and hence makes the theory 
unnecessarily complicated. The astute in doing EFTs is 
in choosing the proper cut-off. Thus with our effective 
Lagrangian in which the lightest degree of freedom inte- 
grated out is the pion, the natural cut-off scale is set by 
the pion mass. We find that the optimal cut-off in our 
case is A ~ 200 MeV as one can see from the results in 
Table 1 and Figures 1 and 2. 

We shall do the calculation to the next-to-leading order 
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(NLO). The potential of the EFT is local and hence of 
zero range in coordinate space, requiring regularization. 
In order to do the calculation algebraically, we choose the 
following form of regularization appropriate to a separa- 
ble potential given by the local Lagrangian: 



(p'|f|p) =5a(p')V^(p'-p)^a(p) 



(1) 



where S'a(p) = ^'(a?) the regulator which suppresses 
the contributions from |p| > A, hmx^o S{x) — 1 and 
limj;3>i S{x) = 0, and V^(q) is a finite-order polynomial 
in q. Up to the NLO, the most general form of V^(q) is 

A-TT .... 

^^(q) = (Co + iC2S'' + D2a^nq'q') , (2) 

where M is the nucleon mass and cr*-' is the rank-two 
tensor that is effective only in the spin-triplet channel, 




(3) 



Note that the coefficients Co, 2 are (spin) channel- 
dependent, and that D2 is effective only in spin-triplet 
channel. Thus we have five parameters; two in ^Sq and 
three in ^Si channel. In principle, these parameters are 
calculable from a fundamental Lagrangian (i.e., QCD) 
but nobody knows how to do this. So in the spirit of 
EFTs, we shall fix them from experiments. Since the ex- 
plicit form of the regulator should not matter [|2|, we 
shall choose the Gaussian form, 



^a(p) =exp ( 



(4) 



where A is the cut-off. (This form of the cut-off func- 
tions, strictly speaking, upsets the chiral counting on 
which we will have more to say later.) The Lippman- 
Schwinger (LS) equation for the wavefunction 1-0), = 

+&V\il}) where \f) is the free wavefunction and 
is the free two-nucleon propagator depending on the to- 
tal energy E, (p'|G°|p 

fmiction (for the potential @)) of the form 
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Ce = OA ('l + ]^aKrKME\ + {VZD2MEfrE, (9) 



with 

GA^Z [Co + (C| + Ss.iDDh)] , 
97 

rA ^ — [2C2 - (C| - ds,iDl)l2] 
"a 

where /„ {n = 2, 4) are defined by 

^n+l foo 



n+1 foo 



(10) 
(11) 



(12) 



With the regulator (Q) , the integrals come out to be I2 = 
-^A3 and h = -Ti-A^ 

The phase shifts can be calculated by looking at the 
large-r behavior of the wavefunction. To do this, it is 
convenient to separate the pole contributions of the inte- 
grals Eqs.(| 0) as 
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which define the functions Ia{E) = A/(- 
H{Ar, ^), both of which are real. Note that H{0, e) = 
1 which makes rA(0) finite, and that lima;3>i H{x, e) = 0. 
The phase shift S takes the form 



p cot 6 
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where the ri{E) is the D/S ratio to be given below (see 
(^)), which vanishes for the ^Sq channel. 

In order to fix the two coefficients Co,2, we compare 
(|l5| ) to the effective-range expansion 
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(16) 

(17) 
(18) 



These are essentially the "renormalization conditions" in 
the standard renormalization procedure. Two important 
observations to make here: (a) We note that there is 
an upper bound of A, Ambx, if one requires that Z be 
positive and that C2 be real. That is, for A > AMax, the 
potential of the EFT becomes non-Hermitian. With a — 
—23.732 fm and — 2.697 fm for the ^Sq channel taken 
from the Argonne Vis potential |l^ (which we take to be 
"experimental"), we find that AMax — 348.0 MeV; (b) the 
value Az=i defined such that Z = 1 when A = Az=i — 
172.2 MeV is quite special. At this point, we have ta = 
and C2 = 0, that is, the NLO contribution is identically 
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zero. This corresponds to the leading-order calculation 
with the A chosen to fit the experimental value of the 
effective range Tg. A similar observation was made by 
Beane et al |^ using a square-well potential in coordinate 
space with a radius R, with R^^ playing the role of A. 

The resulting phase shift with A = Az=i is plotted 
in Fig. 1. We see that the agreement with the result 
taken from the Argonne vis potential is perfect up 
to p ~ 70 MeV. Beyond that, we should expect correc- 
tions from the next-to-next-order and higher-order terms. 
In Fig. 2, we show how the phase-shift for a fixed center- 
of-mass momentum, p — 68.5 MeV varies as the cut-off 
is changed. The solid curve is our NLO result, the dot- 
ted one the LO result (with C2 — 0), and the horizon- 
tal dashed line the result taken from the vis potential 
("experimental"). We confirm that our NLO result is 
remarkably insensitive to the value of A for A > m^. 

It is instructive to compare our result ( |l5| ) with that 
obtained with the dimensional regularization 

1 



PC0t,5|^,„. = ^a-\l + -ar.ME)- 



(19) 



Expanding p cot (5 of (|T^) in ME^ we find that the coeffi- 



cient of the rt-th order term is order of 



This 



increases rapidly with n when a is large, disagreeing 
strongly with the fact that the low-energy scattering is 
well described by just two terms of the effective range 
expansion in (p^). This observation led the authors of |^] 
to conclude that the critical momentum scale at which 
the EFT expansion breaks down is very small for a very 
large a: 



(20) 



We arrive at a different conclusion. With the cut-off 
regularization, the scattering length a is replaced by an 
effective one, a^, that is order of for large a. This 
agrees with the findings of Beane et al |^ and Lepage 
Counting to be order of A~^, the n-th order co- 
efficient now is A^~^", as one would expect on a general 
ground. 

The next quantity to consider is the transition M 1 am- 
plitude for n+p — > d-l-7 and the deuteron structure. For 
the np capture, we need both the ^5*0 scattering wave- 
function and the deuteron wave function. The initial 
state wave function can be written as 
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where p — yME is the center-of-mass momentum and 
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FIG. 1. np ^So phase shift (degrees) vs. the center-of-mass 
(CM) momentum p. Our theory with A = Az=i ~ 172 MeV 
is given by the solid line, and the results from the Argonne 
Vis potential ||10| ("experiments") by the solid dots. 
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FIG. 2. np ^So phase shift (degrees) vs. the cut-ofi A for 
a fixed CM momentum p — 68.5 MeV. The solid curve repre- 
sents the NLO result, the dotted curve the LO result and the 
horizontal dashed line the result from the vis potential ffOl. 



As for the ^Si coupled channel relevant for the fi- 
nal state of the np capture, we use the eigenphase 
parametrization p3) with the ri{E) given by rjlE) = 
— tanei where ei is the mixing angle, 
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The D2 so far undetermined can be fixed by the deuteron 
£1/5' ratio rjd ~ 0.025 ^ at E ^ -Bd with Bd the 
binding energy of the deuteron, 



1 - r;^ -MBd 
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(26) 



Given Co, C'2 and D2 for a given A, all other quantities 
are predictions. The binding energy of the deuteron is 
determined by the pole position, 



iS'{^) + Ia{-i')^ 



1 



GA (1 



^aAVAl^) 



(27) 



with 7 = V MBd- The renormalization procedure is the 
same as for the ^Sq channel. The only difference is that 
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TABLE I. Deuteron properties and the Ml transition am- 
plitude entering into the np capture for various values of A. 



A (MeV) 150 198.8 216.1 250 Exp. [L0| t^is llC] 



Ba (MeV) 
As (fm-*) 
Td (fm) 
Qd (fm^) 
Pd (%) 
fJ-d 

MiB (fm) 



1.799 
0.869 
1.951 
0.231 
2.11 
0.868 
4.06 



2.114 

0.877 

1.960 

0.277 

4.61 

0.854 

4.01 



2.211 
0.878 
1.963 
0.288 
5.89 
0.846 
3.99 



2.389 
0.878 
1.969 
0.305 
9.09 
0.828 
3.96 



2.225 



0, 

1.966(7) 
0.286 

0.8574 



2.225 
0.885 
1.967 
0.270 
5.76 
0.847 
3.98 



As mentioned, the Gaussian cut-ofF brings in terms 
higher order than NLO which to be consistent, would 
require corresponding "counter terms" in the potential 
although our results indicate that the latter cannot be 
significant. The next task is to incorporate pions into 
the picture and go up in energy. This would enable us 
to explore the interplay between the breakdown of EFT 
and the emergence of a "new physics" , an important and 
generic issue currently relevant in particle physics where 
going beyond Standard Model is the Holy Grail. These 
issues will be addressed in a forthcoming publication. 



the value of Az=i that makes Z = 1 does not coincide 
with ArA=o that makes ta = 0. Using a = 5.419 fm 
and re = 1.753 fm for the '^Si channel, we find that 
AMax = 304.0 MeV, Az=i = 198.8 MeV and A^^=o 
216.1 MeV. The resulting (5'-wave and Z?-wave) radial 
wavefunctions of the deuteron are 
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(28) 
(29) 



We now have all the machinery to calculate the 
deuteron properties: the wavefunction normalization fac- 
tor As, the radius r^, the quadrupolc moment Qd and 
the Z?-state probability Pd- The magnetic moment of 
the deuteron fid is related to the Pd through 



fJ-d = fJ-s - ^ { fJ-S 
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Pd 



(30) 



where fJ,s — 0.8798 is the isoscalar nucleon magnetic mo- 
ment. Finally the one-body isovector Ml transition am- 
plitude relevant for n + p ^ d -I- 7 at threshold jij is 



M- 



IB 



dr u{r)uo{r). 



(31) 



The (parameter-free) numerical results are listed in 
Table 1 for various values of the cut-off A. We see 
that the agreement with the experiments (particularly 
for A = 216.1 MeV) is excellent with very little depen- 
dence on the precise value of A. It may be coincidental 
but highly remarkable that even the quadrupole moment 
which as the authors of stressed, the vis potential 
fails to reproduce, comes out correctly. 

We believe to have demonstrated the power of EFTs in 
low-energy nuclear physics, allowing us to be as close as 
one can hope to the fundamental theory in the sense put 
forward in Rcf. [|l]j2|. In particular, it is satisfying that 
the classic np capture process can be completely under- 
stood from a "first-principle" approach. Here the cut-off 
regularization was found to be highly efhcient: With the 
dimensional regularization the Ml matrix element was 
found to be in total disagreement with the result of the 
Argonne vis potential. 
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